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Gluon correlations in the glasma 
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The physics of the initial conditions of heavy ion collisions is dominated by the nonlinear 
gluonic interactions of QCD. These lead to the concepts of parton saturation and the Color 
Glass Condensate (CGC). We discuss recent progress in calculating multi-gluon correlations 
in this framework, prompted by the observation that these correlations are in fact easier to 
compute in a dense system (nucleus- nucleus) than a dilute one (proton- proton). 

o 

CN ■ §1. Introduction 

> 

O ' Bulk particle production in relativistic collisions around midrapidity originates 

. from small x degrees of freedom, predominantly gluons, in the wavefunctions of 

^f) I the colliding hadrons or nuclei. Because of the Inl/x enhancement of soft gluon 

bremsstrahlung this is a dense gluonic system. When the occupation numbers of 
r~| . gluonic states in the wavefunction become large enough, of the order of 1/as, the 

Oh| nonlinear interaction part of the Yang-Mills Lagrangian becomes of the same order 

pi 1' of magnitude as the free part. This leads to the concept of a transverse momentum 

^ . scale Qs, the saturation scale, below which the system is dominated by nonlinear 

interactions. When the collision energy is high enough [x small enough), ^ ^qcd 
and the coupling is weak: we are faced with a nonperturbative strongly interacting 
system with a weak coupling constant. On the other hand, the large occupation 
numbers mean that the system should behave as a classical field. This suggests 
■ a way of organizing calculations that differs from traditional perturbation theory. 

. Instead of developing as a series of powers in gA^ one wants to calculate the classical 

background field A^^ and loop corrections (which are suppressed by powers of g) to 
all orders in gA'^^ . The classical gluon field will then be radiated by the large x 
degrees of freedom, which are treated as effective classical color charges. They can 
be described as random color charges drawn from a classical probability distribution 
Wy [p] that depends on the rapidity cutoff y = Inl/x separating the large and small 
^ . X degrees of freedom. The dependence of Wy[p] on y is described by a Wilsonian 

^ [ renormalization group equation known by the acronym JIMWLK. This picture of 

the high energy wavefunction is referred to as the Color Glass condensate (CGC, for 
reviews see e.g.^). 

The role of Wy [p] is analoguous to the conventional parton distribution function; 
it is a nonperturbative quantity whose dependence on one of the kinematical variables 
of the process is described by a weak coupling renormalization group equation. In the 
case of pdf 's the appropriate degrees of freedom are individual partons with a definite 
momentum, whereas in the case of the CGC they are color charges resulting from 
interactions of many partons. The distributions are, like pdf's, not (complex) 
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Fig. 1. Left: Building block, Lipatov vertex coupled to two sources. Right: combinatorics of the 
sources. The combinatorial problem is to connect the dots on the upper and lower side (left- 
and right moving sources) pairwise. 



wavefunctions but real and can be interpreted as probability distributions. This is 
guaranteed by factorization theorems^ at large x and the process one is studying. 
Factorization can be understood as a statement that one has found the right set 
of degrees of freedom in which one can compute physical observables from only the 
diagonal elements of the density matrix of the incoming nuclei. 

The term glasmd^ refers to the coherent, classical field configuration resulting 
from the collision of two objects described in the CGC framework. The glasma fields 
are initially longitudinal, whence the "glasma fiux tube"'^''3' picture. More impor- 
tantly for computing multigluon correlations, they are boost invariant (to leading 
order in the QCD coupling) and depend on the transverse coordinate with a charac- 
teristic correlation length 1/Qs- There are several signals in the RHIC datd^'l^ that 
point to strong correlations originating from the initial stage of the collision. The 
glasma fields provide a natural framework for understanding these effects, although 
much work is still left to do in understanding the interplay with purely geometrical 
effects from the fluctuating positions of the nucleons in the colliding nuclei.'®' 

We shall flrst make some general observations on computing multigluon correla- 
tions in the glasma, arguing in Sec. [2] that they are in some sense simpler to compute 
in a collision of two dense, saturated nuclear wavefunctions than in the dilute limit. 
We shall then, in Sec. [3] discuss one application of these ideas to computing the mul- 
tiplicity distribution of gluons in the collision before moving to the leading Inl/x 
rapidity dependence of the correlation in Sec. [H The application of these ideas to 
understanding correlations observed in the experimental data are discussed in Sec. [5l 

§2. Multigluon correlations in the glasma 

The gluon fields in the glasma are nonperturbatively strong, ^ 1/g. This 
means that the gluon multiplicity is ~ l/os- For a fixed configuration of the 
classical color sources it is well known that the multiplicity distribution of produced 
gluons is Poissonian, i.e. {N'^) — (-/V)^ = {N). In this case the correlations and fluc- 
tuations in the gluon multiplicity are all quantum effects that appear only starting 
from the one-loop level, i.e. suppressed by a power of the coupling constant a^- The 
computation in the CGC framework does not end here, however. To calculate the 
moments of the gluon multiplicity distribution one must first calculate the gluon 
spectra for fixed configuration of the color charges p and then average over the prob- 
ability distribution Wy[/9(a;^)]. For the nth moment of the multiplicity distribution, 
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i.e. an n-gluon correlation, the leading order result is 



( 
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(2-1) 



[p] 



where the subsrcipt "LO" refers to the single gluon spectrum evaluated from the 
classical field configuration corresponding to a fixed configuration of color charges. 
This averaging, even after the subsequent subtraction of the appropriate discon- 
nected contributions, introduces a correlation already at the leading order in as, i.e. 
enhanced by an additional l/og compared to the quantum correlations. A natural 
example is the negative binomial distribution that we shall discuss below, whose 
variance is {N'^) — {N)'^ = {N)'^/k + (N). One must emphasize here that although 
these contributions arise as formally classical correlations in the effective theory that 
is the CGC, they are physically also quantum effects, where the weak coupling is 
compensated by a large logarithm of the energy that has been resummed into the 
probability distribution Wj^[/9(a;x)]. In this sense the leading correlations are present 
already in the wavefunctions of the colliding objects. 



We can then apply this formalism to the calculation of the probability distri- 
bution of the number of gluons in the glasma.^' We shall assume the "AA" power 
counting of sources p that are parametrically strong in g, but nevertheless work to 
the lowest nontrivial order in p. Formally this would correspond to a power counting 
p ~ with a small e > 0. In this limit, as we have discussed, the dominant contri- 
butions to multiparticle correlations come from diagrams that are disconnected for 
fixed sources and become connected only after averaging over the color charge config- 
urations. The corresponding two gluon correlation function was computed in Ref.® 
and generalized to a three gluons in Ref.'^ We shall here sketch the derivation'^ of 
the general n- gluon correlation in this simplified limit. 

Working with the MV model Gaussian probability distribution 



computing the correlations and the multiplicity distribution in the linearized ap- 
proximation is a simple combinatorial problem. Each gluon is produced from two 
Lipatov vertices (see fig. [1] left), one in the amplitude and the other in the complex 
conjugate. The combinatorial factor is obtained by counting the different ways of 
contracting the sources pairwise (see fig. [1] right). The dominant contributions are 
the ones that have, in the dilute limit, the strongest infrared divergence which is 
regulated by the transverse correlation scale of the problem, Qg. When integrated 
over the momenta of the produced gluons one obtains the factorial moments of the 
multiplicity, which define the whole probability distribution. It can be expressed 
in terms of two parameters, the mean multiplicity n, and a parameter k describing 



§3. Multiplicity distribution 
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Fig. 2. Numerical evaluation of the two-gluon correlation strength in the MV model as a function 
of the two transverse momenta p± and q±. Left: the correlation strength on the "near side" 
Ifp — 'fiq\ < tt/S, right: the "away side" — ipq\ > n/2. 



the width of the distribution. The resuh of the combinatorial exercise is that the 
qth factorial moment niq (defined as (A^"^) minus the corresponding disconnected 
contributions) is 



iriq = {q — l)\k 



with k 



2tt 



and n 



1 



a. 



(3-2) 



Here Sx_ is the transverse area of the system. These moments define a negative 
binomial distribution (NBD) with parameters k and n. The NBD has been known as 
a purely phenomenological observation in high energy hadron and nuclear collisions 
already for a long time. Also the numerical magnitude of the parameter k obtained 
from the saturation scale agrees very well with both pp and AA data.'^''^ 

In terms of the glasma flux tube picture this result has a natural interpretation. 
The transverse area of a typical flux tube is l/Qg, and thus there are QsS± = Np^ 
independent ones. Each of these radiates particles independently into Nc'^ — 1 color 
states in a Bose-Einstein distribution (see e.g.^^). A sum of ~ A^ft(-^c^ — 1) in- 
dependent Bose-Einstein-distributions is precisely equivalent to a negative binomial 
distribution with parameter k. 

This calculation predicts that the k parameter should increase with energy as 
^ ~ Qs ~ y/s^- This is ideed true in the heavy ion data; k is reported to increase 
from = -462 GeV to ^200 GeV in by the PHENIX cohaboration.i^J However, 
the interpretation of the heavy ion data is complicated by the purely geometrical 
fluctuations from the different impact parameters probed in one centrality bin. 

The proton-proton collision system is much smaller and fluctuations at lower 
energies are still mostly dominated by the dilute edge of the collision system, which 
has a Poissonian nature (i.e. k — )• oo). Our calculation formally assumes QsS± ^ 1, 
and we expect the growing behavior of k with energy to eventually take over at 
high enough energy. Since the fluctuations are more dominated by the edge region 
than the mean multiplicity, it is natural to expect this transition to a genuine high 
energy regime to be visible later in the fluctuations than it is in the average. At 
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Fig. 3. The parameter k in pp collisions measured by the UA5p' and ALICe!^ For UA5 the 
NBD fit is done by the experiment. The fit to ALICE data is done by us and does not take into 
account the correlated errors in the experimental data points. 

LHC energies (see the data in Ref.'^) the transition to an increasing k seems to be 
starting as shown in fig. O 



§4. Rapidity dependence 

The general discussion of Sec.[2]on the different nature of multigluon correlations 
in the "AA" case applies also to the rapidity dependence. Until now we have only 
been discussing gluon production in a rapidity interval smaller than l/ag. For this 
we needed only the correlations between the color charges p{x±) measured at this 
same rapidity. To understand the rapidity dependence of the correlations one needs 
also the correlation between color charges at different rapidities, {pyix±)py'{y^)). 
Also this information is contained in the JIMWLK renormalization group evolution, 
at least to leading Inl/x accuracy.'^ An intuitive description of the resulting cor- 
relations is provided by the formulation of JIMWLK as a Langevin equation in the 
space of Wilson lines formed from the color charges. In this picture the evolution 
proceeds in individual trajectories along an increasing rapidity. 

A first attempt of a realistic estimate of the rapidity dependence of two-gluon 
correlations is performed in Ref.'^ Evaluating the two gluon correlation in a dilute 
limit in a /c^-factorized approximation, but keeping the general structure resulting 
from the JIMWLK evolution leads to the following expression: 



as 



+ ^A2iyq^k^)'^AAyp^p± - k±)\'pAAyq^Q± + k±) + '^AAyq^Q± - k±)\ }■ (4-i) 



Note the very different structure of this correlation compared to one where the glu- 
ons would be produced from the same diagram for fixed sources. The two gluon 
correlation function is proportional to the product of four unintegrated gluon distri- 
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Fig. 4. Left: Comparison of a two-particle correlation computed using eq. (|4-l|l supplemented with 
a short-range correlation contribution from PYTHIA with PHOBOS data. Right: Rapidity 
correlation at LHC energies fcx-factorization approximation. 



butions, with three of them evaluated at the rapidity of one of the produced gluons 
and only one at the other. This structure is a direct consequence of the nature of 
JIMWLK evolution. 

§5. The ridge in nuclear and proton collisions 

The "ridge" is a feature on the "near side" of the two particle correlation, around 
A(j) ~ azimuthal separation between the two particles. It extends up to high val- 
ues of Ar], the pseudorapidity difference between the two particles. This feature 
was first seen in nucleus-nucleus collisions at RHICl^^^ The ridge was seen in 
high multiplicity (central) events in nucleus-nucleus collisions in Cu-|-Cu collisions 
at ^ = 62.4 GeV and in Au-hAu collisions at ys = 200 GeV. The STAR detector 
observed this correlation for both p_L-triggerecP^ and untriggerecF^ pair correla- 
tions in the whole STAR TPC acceptance of Ar] < 2. The PHOBOS experimenlPS 
observed the p_L-triggered correlation at much larger rapidity separations, initially 
up to At] ~ 4, extended more recentl}^^ to Ar] ~ 5. The long range correlation 
structure disappears for lower multiplicity peripheral events in nucleus-nucleus colli- 
sions and are also absent in deuteron-gold and proton-proton "control" experiments 
at the same energies. 

The "standard" glasma explanationP of the azimuthal structure of the ridge 
relies on a collimation effect from radial flow combined with the long range rapidity 
correlation from the boost invariant color fields. The correlation computed from 
eq. (|4-1|) is compared to PHOBOS data in fig. HI The A:^-factorized approximation 
gives a very inaccurate description of the gluon spectrum in the transverse momen- 
tum regime p± ~ Qs where the bulk of the particles are produced .'^S' Equation ()4-ip 
has also been derived in the approximation, true only in the linearized case, that 
the unequal rapidity correlation of two color charge densities is equal to the unin- 
tegrated gluon distribution at the smaller one of these rapidities. As of yet there 
is no calculation of how much this approximation is violated in the full JIMWLK 
evolution (see, however ,123). The results presented in fig. |1] are therefore not the 
final word on the subject, alhough it is reassuring that such a simple approximation 
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Fig. 5. Two-gluon correlation strength in the MV model. Both figures show the same data; the 
large delta function peaks have been removed from the one on the right to illustrate the structure 
in the rest of phase space. 



is shown in Fig. [2j It confirms the expectations of Ref.^' on the magnitude of the 
correlation. In contrast to the approximation made in Ref.^^ there is, however, 
some dependence on the momenta which is crucial for the discussion of the 

correlation in smaller systems where transverse flow is smaller or absent. 

More recently the CMS collaboration has reported the observation of a "ridge" 
structure also in proton-proton collisions. The ridge is seen only in the moderate 
P± ) Q± range and systematically vanishes for p±,q± < 1 GeV and p±,q± ^ 3 GeV. 
Although the explanation in terms of the collimating effect of the transverse flow is 
most likely present to some degree in high multiplicity proton-proton collisions at 
LHC energies, this observation points to a ridge-like effect being present already in 
the initial scattering process. This is indeed the picture resulting from the correla- 
tions in the glasma. 

In terms of the fc^-factorized assumption (|4Tp the CMS data was recently dis- 
cussed in Ref.l^ The dependence of the CMS ridge on transverse momentum was 
found to qualitatively agree with the CMS result. In the full nonperturbative calcu- 
lation this same collimation effect already in the initial stage is visible in fig. [5l where 
the correlation is plotted as a function of — q_|_| and \pj_ + q±\- The production 
of gluons from a coherent classical field results in an enhanced near side correlation 
for momenta and q that are close to parallel, even in the absence of transverse 
flow. 



seems to agree rather well with the experimental result. A numerical evaluatiorJ^ 
of the second moment of the distribution, parametrized in terms of 





§6. Conclusion 



Most experimental observables do not probe the glasma initial state of directly, 
because the system goes through a complicated time evolution before the hadroniza- 
tion stage. A good candidate for an experimental probe giving direct access to the 
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initial state is provided by different kinds of correlation measurements. These have 
indeed been a focus of both experimental and theoretical activity recently. We have 
argued here that the glasma picture of the initial stages of a heavy ion collision is 
the natural framework to understand the origin of these correlations. 
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